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The square lattice is perhaps the simplest cellular structure. In this work, however, we investigate 
the various structural and topological properties of the kinetic and stochastic counterpart of the 
square lattice and termed them as kinetic square lattice (KSL) and weighted planar stochastic 
fvj ' lattice (WPSL) respectively. We find that WPSL evolves following several non-trivial conservation 
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5_j ■ laws, x™~ 1 y™ = const. V n, where x, and ?/; are the length and width of the ith block. 

f> |' The KSL, on the other hand, evolves following only one conservation law, namely the total area, 

although one find three apparently different conserved integrals which effectively the total area. We 
show that one of the conserved quantity of the WPSL obtained either by setting n — 1 or n = 4 can 
be used to perform multifractal analysis. For instance, we show that if the ith block is populated 
with either pt ~ a;f or pi ~ yf then the resulting distribution in the WPSL exhibits multifractality. 
Furthermore, we show that the dual of the WPSL, obtained by replacing each block with a node 
' pH ' ' at its center and common border between blocks with an edge joining the two vertices, emerges 

as a scale-free network since its degree distribution exhibits power-law P(k) ~ fc~ 7 with exponent 
. 7 = 5.66. It implies that the coordination number distribution of the WPSL is scale-free in character 

as we find that P(k) also describes the fraction of blocks having k neighbours. 

! PACS numbers: 89.75.Fb,02.10.Ox,89.20.Hh,02.10.Ox 

£ 

The formation of planar cellular structures have attracted considerable interest among scientists in general and 
physicists in particular because cellular structures are ubiquitous in nature. Especially, the class of non-equilibrium 
£h ■ systems with apparent disordered patterns have come under close scrutiny. Examples of cellular structures include 
acicular texture in martensite growth, tessellated pavement on ocean shores, agricultural land division according to 
ownership, grain texture in polycrystals, cell texture in biology, soap froths and so on [J-Q. Most existing theoretical 
models which can mimic such systems are typically formed by tessellation, tiling, or subdivision of a plane into domains 
of contiguous and nonoverlapping cells. For instance, Voronoi lattice is formed by partitioning a plane into mutually 
exclusive convex polygons, Apollonian packing on the other hand is formed by tiling a plane into contiguous and 
non-overlapping disks, etc. These models have found widespread applications in physics and biology [H, H|. Significant 
progress has already been made in understanding various properties of both Voronoi lattice and Apollonian packing. 

A fact worth mentioning is that most of the existing models for generating cellular structures have one thing in 
common and that is: the number of neighbours of a cell never exceeds the number of sides it has. Another important 
property of the existing cellular structures is that they have small mean coordination number where it is almost 
impossible to find cells that have significantly higher or fewer neighbours than the average number of neighbours. 
However, geologists and soil scientists observe an abundance of cellular patterns ranging from topographical ground 
formations and underwater stone arrangements where number of neighbours may exceeds the number of sides of the 
cell. Typically, these structures emerge through evolution and in the process each cell may acquire more neighbours 
than the number of sides the corresponding cell has. Such disordered structures can be of great interest in physics 
provided they have some global topological and geometrical properties since they can mimic kinetic disordered medium 
on which one can study percolation and random walk problems. 

Perhaps, the square lattice is the simplest example of the cellular structure where every cell has the same size and 
the same coordination number. Its construction starts with an initiator, say a square of unit area, and a generator that 
divides it into four equal parts. In the next step and steps thereafter the generator is applied to all the available blocks 
which eventually generates a square lattice. In this talk, we intend to address the following questions. Firstly, what if 
the generator is applied to only one of the available block at each step by picking it preferentially with respect to the 
areas? Secondly, what if we use a modified generator that divides the initiator randomly into four blocks instead of 
four equal parts and apply it to only one of the available block at each step, which are again picked preferentially with 
respect to their respective areas Q? Our primary focus, however, will be on the later case which results in the tiling 
of the initiator into increasingly smaller mutually exclusive rectangular blocks. The process is so simple that even a 
few steps by hand on a piece of paper can lead to the conclusion that the blocks in the resulting lattice may have more 
neighbours than the number of sides of a cell. We term the resulting structure as weighted planar stochastic lattice 



FIG. 1: A snapshot of the weighted planar stochastic lattice containing 30001 blocks. 



(WPSL) since the spatial randomness is incorporated by the modified generator and also the time is incorporated in 
it by the sequential application of the modified generator. The definition of the model may appear too simple but 
the results it offers are found to be far from simple. To illustrate the type of systems expected, we show a snapshot 
of the resulting weighted planar stochastic lattice taken during the evolution (see figure 1). We intend to investigate 
its topological and geometrical properties in an attempt to find some order in this seemingly disordered lattice. 

The model in question can also describe the kinetics of fragmentation in two dimensions as it can be defined as 
follows At each time step one may assume that a seed is being nucleated randomly on the initiator and hence the 

greater the area of the block, the higher is the probability that it contains the seed. Upon nucleation two orthogonal 
cracks parallel to the sides of the block are grown until intercepted by existing cracks which divides it randomly into 
four rectangular fragments. In reality though, fragments are produced in fracture of solid objects by propagation of 
interacting cracks resulting in the fragments of arbitrary size and shape which makes it a formidable mathematical 
problem. However, the present model can be considered as the minimum model which should be capable of capturing 
the essential generic aspects of the underlying mechanism. In addition, Ben-Nairn and Krapivsky have noted that 
the WPSL can also describe kinetics of martensite formation. A clear testament to this can be found if one compares 
figure 1 of our work and figure 2 of the work of Rao et al. in the context of martensite formation [T], [T(| • Yet another 
application, albeit a little exotic, is the random search tree problem in computer science [Til ]. 

A truly disordered lattice where the coordination number disorder and the block size disorder are introduced 
naturally is perhaps long over due. On the other hand, where there is a disorder physicists have the natural tendency 
to look for an order in it. To this end, we invoke the idea of multifractality to quantify the size disorder and the idea 
of scale-invariance to quantify the coordination number disorder. Firstly, we characterize the blocks of the lattice 
by the size of their respective length x and width y. Furthermore, we identify each blocks by labelling them as 
i = 1,2, ...,N and the corresponding length and width by (x\,y\), (2:2,2/2), • (xn,Hn) respectively. We then show 
that the dynamics of the WPSL is governed by infinitely many conservation laws, namely the numerical value of 
x i ~ 1 yt^ n 1 remains the same regardless of the lattice size N and the value of n where n = 2 corresponds to 
the conservation of total area. Of the non-trivial conservation laws, we single out either YJ^ yf by setting n = 1 or 
x| by setting n — 4 as we can use it to perform multifractal analysis. For instance, we show that if the blocks 
are populated according to the cubic power of their own length (or width) then the distribution of the population in 
the WPSL emerges as multifractal. On the other hand, if the blocks are characterized by the coordination number 
or the number of neighbours k with whom the block have common sides then we show that the coordination number 
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distribution function p(k), the probability that a block picked at random have coordination number k, have power-law 
tail. A power-law distribution function is regarded as scale-free since it looks the same regardless of the scale we use 
to look at it. 

The organization of this talk is as follows. In section 2, we give the exact algorithm of the model. In section 3, 
we give the geometric properties of the WPSL in an attempt to quantify the annealed size disorder and at the same 
time we proposed kinetic square lattice in an attempt to look for a possible origin of multifractality. In section 4, we 
discuss the various structural topological properties of the WPSL and its dual are discussed in order to quantify the 
annealed coordination number disorder. Finally, section 5 gives a short summary of our results. 
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STEP 1 



STEP 5 & 6 (dotted lines) 



FIG. 2: Schematic illustration of the first few steps of the algorithm. 



II. ALGORITHM OF THE WEIGHTED PLANAR STOCHASTIC LATTICE (WPSL) 

Perhaps an exact algorithm can provide a better description of the model than the mere definition. In step one, the 
generator divides the initiator, say a square of unit area, randomly into four smaller blocks. The four newly created 
blocks are then labelled by their respective areas oi, (12,(13 and 04 in a clockwise fashion starting from the upper left 
block (see Fig. 2). In step two and thereafter only one block is picked at each step with probability equal to their 
respective area and then it is divided randomly into four blocks. In general, the jth step of the algorithm can be 
described as follows. 

(i) Subdivide the interval [0, 1] into (3j — 2) subintervals of size [0, ai], [01, a\ + 02], Ei£7 3 a »' 1] eacn of which 
represent the blocks labelled by their areas a\,ai, ...,a(3j_2) respectively. 

(ii) Generate a random number R from the interval [0, 1] and find which of the (3i — 2) sub-interval contains this 
R. The corresponding block it represents, say the pth block of area a p , is picked. 

(iii) Calculate the length x p and the width y p of this block and keep note of the coordinate of the lower-left corner 
of the pth block, say it is {xi ow ,yi ow ). 

(iv) Generate two random numbers xr and yji from [0, x p ] and [0, y p ] respectively and hence the point (xr+xi ow , yn+ 
yiow) mimicking a random nucleation of a seed in the block p. 

(v) Draw two perpendicular lines through the point (xr + xi ow ,yn + yi ow ) parallel to the sides of the pth block 
mimicking orthogonal cracks parallel to the sides of the blocks which stops growing upon touching existing 
cracks and divide it into four smaller blocks. The label a p is now redundant and hence it can be reused. 

(vi) Label the four newly created blocks according to their areas a p , a(3j_i), 03-,- and a/sj+u respectively in a 
clockwise fashion starting from the upper left corner. 

(vii) Increase time by one unit and repeat the steps (i) - (vii) ad infinitum. 
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III. GEOMETRIC PROPERTIES OF WPSL 



In general, the distribution function C(x, y; t) describing the blocks of the lattice by their length x and width y 
evolves according to the following kinetic equation || [t| 

= -C(x,y;t) / dxidyiF{xi,x-xi,yi,y-yi)+ (1) 



dt 



o Jo 



oo fOO 



4/ / C(x 1 ,y 1 ;t)F(x,x 1 - x,y,yi - y)dx 1 dy 1: 

J x J y 

where kernel F(xi, x 2 , yi, y 2 ) determines the rules and the rate at which the block of sides (x\ + X2) and (yi + y 2 ) is 
divided into four smaller blocks whose sides are the arguments of the kernel. The first term on the right hand side of 
equation (fTJ) represents the loss of blocks of sides x and y due to nucleation of seed of crack on one such block from 
which mutually perpendicular cracks are grown to divide it into four smaller blocks. Similarly, the second term on the 
right hand side represents the gain of blocks of sides x and y due to nucleation of seed of crack on a block of sides x\ 
and yi ensuring that one of the four new blocks have sides x and y. Let us now consider the case where the generator 
divides the initiator randomly into four smaller rectangles and apply it to only one of the available squares thereafter 
by picking preferentially with respect to their areas. It effectively describes the random sequential nucleation of seeds 
with uniform probability on the initiator. Within the rate equation approach this can be ensured if one chooses the 
following kernel 

F(x 1 ,x 2 ,yi,y 2 ) = 1 (2) 

Substituting it into equation (|T|) we obtain 

dC(x,y;t) 



x Jy 



- -xyC(x,y;t)+4 / / C(xi,yi;t)dxidyi. (3) 



The coefficient xy of C(x, y; t) of the loss term implies that seeds of cracks are nucleated on the blocks preferentially 
with respect to their areas which is consistent with the definition of our model. 
Incorporating the 2-tuple Mellin transform given by 
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m — l„,n— 1/ 



M(m,n;t)= / x m - l y n - l C{x,y-t)dxd yi (4) 
Jo Jo 



in equation Q we get 

dM(m,n;t) / 4 



dt 



( lW(m + l,n + l;i). (5) 



Iterating equation ([5]) to get all the derivatives of M(m, n; t) and then substituting them into the Taylor series 
expansion of M(m,n;t) about t = one can immediately write its solution in terms of generalized hypergeometric 
function [l2[ 

M(m,n;t) = 2-p2fa_|_,a_;m, n\— tj, (6) 
where M(m, n; t) — M(n, m; t) for symmetry reason and 

-/m — n\ 2 13 



m + n 



a± = ± 



One can see that (i) M (1, 1; t) = 1 + 3t is the total number of blocks N(i) and (ii) M(2, 2; t) — 1 is the sum of areas 
of all the blocks which is obviously a conserved quantity. Both properties are again consistent with the definition of 
the WPSL depicted in the algorithm. The behaviour of M(m, n; t) in the long time limit is 

M{m,n;t) - t~ a - . (8) 

Thus, in addition to the conservation of total area the system is also governed by infinitely many non-trivial conser- 
vation laws as it implies 



M(n,A/n;t) ~ constant Vn. 



(9) 
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FIG. 3: The plots of x™ 1 yf^ n 1 vs iV for n = 3, 4, 5 are drawn using data collected from one realization. 
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FIG. 4: The plots of vs iV for four different realizations shows that the numerical value is different at every independent 

realization. 



We used numerical simulation to verify equation or its discrete counterpart x 7~ 1 Vi 1 ^ we l & bel all the 
available blocks as i = 1,2, ...,N. We found that the analytical solution is in perfect agreement with the numerical 
simulation which we performed based on the algorithm for the WPSL model (see figure 3) . 

We now find it interesting to focus on the distribution function n(x,i) = J Q C{x,y 1 t)dy that describes the con- 
centration of blocks which have length x at time t regardless of the size of their widths y. Then the qth moment of 
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n(x,t) is defined as 



M q {t) 



x q n(x, t)dx. 



(10) 



Appreciating the fact that M q {t) = M(q + 1,1; t) and using equation ([5J one can immediately write its solution 



M,(t) ~ t*VTO8-(fl+2)}/ 2i (u) 

Note that for symmetry reason it does not matter whether we consider the qth moments of n(x, t) or that of n(y, t) 
since we have M(q + 1,1; t) = M(l,q + l;t). According to equation (fTTI) the quantity Ms(t) and hence J2i x i 
or £\ yf is a conserved quantity. However, although M 3 (t) remains a constant against time in every independent 
realization the exact numerical value is found to be different in every different realization (see figure 4). It clearly 
indicates lack of self-averaging or wild fluctuation. This is supported by our analytical solution as we find that 



J °° x q n(x, t)dx 



^<x> q = ( 



xn(x,t)dx\Q 



J °° n(x,t)dx \ J °° n(x,t)dx ' ' 

which suggest that a single length scale cannot characterize all the moments of the distribution function n(x,t). 



(12) 




FIG. 5: The plots of r(q) vs q to show its slope as a function of q varies. 

Of all the conservation laws we find that M^(t) — J2i x i ^ s a special one since we can use it as a multifractal 
measure consisting of members pi, the fraction of the total measure pi = #f/X)i a '?) distributed on the geometric 
support WPSL. That is, we assume that the «th block is occupied with cubic power of its own length X{. The 
corresponding "partition function" of multifractal formalism then is (l3j 



i 

Its solution can immediately be obtained from equation (jlll) to give 



Z q (t) - t { v /9q2+16 ~ {3q+2)}/2 



(13) 



(14) 



Using the square root of the mean area S(t) = \J M(2, 2; t)/M(l, 1; t) ~ t~ x / 2 as the yard-stick to express the partition 
function Z q gives the weighted number of squares N(q, S) needed to cover the measure which we find decays following 
power- law 



(15) 
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where the mass exponent 



r(q) = VV + 16- (3g + 2). 
The non-linear nature of r(<?), see figure 5 for instance, suggests that the gap exponent 

A = r(q) - r(q - 1) 



(16) 



(17) 



is different for every q value. It implies that we require an infinite hierarchy of exponents to specify how the moments 
of the probabilities {p}s scales with S. Now if we choose q = then it gives an estimate of the number of squares 
N(0, S = N(6) of sides S we need to cover the support on which the members of the population is distributed. We 
find that N(8) scales as 



N(q = 0,6) ~ 6- Tt >°\ 



(18) 



where r(0) = 2 is the Hausdorff-Besicovitch dimension of the support [14]. On the other hand, if we choose q = 1 
we have Z± — ^Zi? 5 * = const, and hence we must have r(l) = 0. This is indeed the case according to equation (|16|l . 
Therefore, r(0) = 2 (the dimension of the support) and r(l) = (required by the normalization condition) are often 
considered as the first self-consistency check for the multifractal analysis. 

The role of q in the partition function Z q can be best understood by appreciating the fact that the large values 
of q favours contributions in the Z q from the blocks with relatively high values of pi since pf » Pj with pi > pj 
provided q >> 1. On the other hand, q << — 1 favours contributions in the Z q from those blocks which are occupied 
with relatively low values of the measure pi. For further elucidation we find it worthwhile to consider the slope of the 
curve r(q) vs q which is given by 



dq «->oQ>?)]n*- 



Now r(q) vs q curve has the maximum slope in the limit q 

dr(q) _ 
dq 

It implies that the right hand side of equation (|19p is dominated by p„ 
we have 



-oo and hence we can write 



(19) 



(20) 



dr(q) 



dq 



q^ — co 



hm — — — 

8-yo In 5 



the minimum of pi , in the sum and hence 



-awx, (21) 



and hence we get 



A similar argument in the limit q 
by 



Pmxn - S a —. (22) 
oo leads to the conclusion that the minimum slope of the r(q) vs q curve is given 



dr(q) 



dq 



= — lim 

q— > + oo 5— >0 



lnp„ 



In 5 



where p max is the largest value of p which gives the minimum value of a and hence 

Fmax u 

So, in general we can write 

dr(q) 



dq 



= -a(q), 



(23) 



(24) 



(25) 



where the exponent a(q) is known as the Lipschitz-Holder exponent and 



p~5 a . 



(26) 
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FIG. 6: The /(a) spectrum. 



We now perform the Legendre transformation of the mass exponent r(q) by using the Lipschitz-Holder exponent 
given equation (|25|) as an independent variable to obtain the new function 



f(a) = qa + T(q). 
Replacing r(q) in equation (|15[) in favour of f(a) we find that 

N(q(a),S) ~ lim 5 qa ~ f{a) . 

8—>0 



(27) 



(28) 



On the other hand, using p ~ <5 Q in the expression for partition function Z g = £\ p q and replacing the sum by integral 
while indexing the blocks by continuous Lipschitz-Holder exponent a as variable with a weight p(a) we obtain 

N(q{a),S) ~ J p(a)daN(a,S)S qa , (29) 

where N(a,S) is the number of squares of side S needed to cover the measure indexed by a. Comparing Eqs. (|28[) 
and (gU) we find 



N(a,S) ~(T /(Q) . 

It implies that we have a spectrum of spatially intertwined fractal dimensions 



/(«(«)) 



1G 



VV + 16 



(30) 



(31) 



are needed to characterize the measure. That is, the size disorder of the blocks are multifractal in character since the 
measure {p a } is related to size of the blocks. That is, the distribution of {p a } in WPSL can be subdivided into a 
union of fractal subsets each with fractal dimension f(a) < 2 in which the measure p a scales as S a . Note that f(a) 
is always concave in character (see figure 6) with a single maximum at q = which corresponds to the dimension of 
the WPSL with empty blocks. 

On the other hand, we find that the entropy S(S) = — '^2 i Pi lnp^ associated with the partition of the measure on 
the support (WPSL) by using the relation ^2 t pf ~ 8^ T ^ in the definition of S(S). Then a few steps of algebraic 
manipulation reveals that S(S) exhibits scaling 



S(6) = ln<r Q W 



(32) 
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where the exponent a\ = | obtained from 



a(q) = - dT(q)/dq\ q . (33) 

It is interesting to note that ct(l) is related to the generalized dimension D q , also related to the Renyi entropy 
H q (p) = — 3y hi^jpf in the information theory, given by 



D a — lim 



1 InE^fU^ (34) 



1 In S 



which is often used in the multifractal formalism as it can also provide insightful interpretation. For instance, 
Dq = t(0) is the dimension of the support, D\ = a\ is the Renyi information dimension and D 2 is known as the 
correlation dimension [l^, [l|| . Multifractal analysis was initially proposed to treat turbulence but later successfully 
applied in a wide range of exciting field of research. For instance, it has been found recently that the wild fluctuations 
of the wave functions at the Anderson and the quantum Hall transition can be best described by multifractality 
[l7l [l8j . Recently, though it has got renewed momentum as it has been found that the probability density function 
at the Anderson and the quantum Hall transition exhibits multifractality since in the vicinity of the transition point 
fluctuations are wild - a characteristic feature of multifractal behaviour [18( . 

In an attempt to understand the origin of multifractality we now consider the case where the generator divides the 
initiator into four equal blocks instead of randomly into four blocks. If the generator is applied over and over again 
thereafter to only one of the available squares by picking preferentially with respect to their areas then it results in 
the kinetic square lattice (KSL). Within the rate equation approach it can be described by the kernel 

F(xi, £2,2/1,2/2) = {xi + £2X2/1 + V2)5{xi - x 2 )5(y 1 - y 2 ). (35) 
and hence the resulting rate equation can be obtained after substituting it in equation ([T]) to give 

dC(x,y;t) _ 1 xyC ^ x y .^ + 4?xyC(2x,2y;t). (36) 



dt 4 
Incorporating equation (|4]) in equation ([36j) yield 

dM(m,n;t) _ (I 



dt V4 2 m+n 



M(m+l,n + l;t). (37) 



To obtain the solution of this equation in the long-time limit, we assume the following power-law asymptotic behaviour 
of M (to, n; t) and write 

M(m, n; t) - A{m, n)t e< - m+n ^ , (38) 

with 9(4) — since total area obtained by setting m — n = 2 is an obvious conserved quantity. Using it in equation 
(l37l) yields the following difference equation 



e(m + n + 2) = 6(m + n) - 1. (39) 

Iterating it subject to the condition that 0(4) = gives 

M(m,n;t) - t - m+ r i . (40) 

Apparently it appears that in addition to the conservation of the total area M(2,2;i), we find that the integrals 
M(3, and M(l,3;£)) are also conserved. Interestingly, all the three integrals M(2, 2;t), M(3,l;i) and M(l,3;t) 
effectively describe the same physical quantity since all the blocks are square in shape and hence 



N N N 

i=l i=l i=l 



&i =2»i =£)aWi- ( 41 ) 



Therefore, in reality the system obeys only one conservation law - conservation of total area. 

We again look into the qth moment of n(x,t) using equation (fT0|) and appreciating the fact that M q (t) equal to 
M(q + 1, 1; t) or M(l, q + 1; t) we immediately find that 



Af,(t)~t-V. 



(42) 
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Unlike in the previous case where the exponent of the power-law solution of the M q (t) is non- linear, here we have an 
exponent which is linear in q. It immediately implies that in the case kinetic square lattice 

< x q >=< x > 9 , (43) 

and hence a single length-scale is enough to characterize all the moments of n(x,t). That is, the system now exhibits 
simple scaling instead of multiscaling. Like before let us consider that each block is occupied with a fraction of the 
measure equal to square of its own length or area pi = Y^f x t an d hence the corresponding partition function is 

N 

Z q =Y J V\ = M M {t)- (44) 

i 

Using equation f|42p we can immediately write its solution 

Z q (t)~t- 2 -^. (45) 

Expressing it in terms of the square root of the average area S ~ i^ 1 / 2 gives the weighted number of square N(q, 5) 
of side S needed to cover the measure which has the following power-law solution 

N(q,8) ~6- {2 - 2q \ (46) 

where mass exponent r(q) is 

r(q) =2-2q. (47) 
The Legendre transform of the mass exponent is a constant 

f(a) = 2, (48) 

and so is the generalized dimension 

D q = 2. (49) 

We thus find that if the generator divides the initiator into four equal square and we apply it thereafter sequentially 
then the resulting lattice is no longer a multifractal. The reason is that the distribution of the population in the 
resulting support in this case is uniform. The two model therefore provides a unique opportunity to look for possible 
origin of multifractality. The two models discussed in this talk differs only in the definition of the generator. In the case 
when the generator divides the initiator randomly into four blocks and we apply it over and over again sequentially 
then we have multifractality since the underlying mechanism in this case is governed by random multiplicative process. 
This is not, however, the case if the generator divides the initiator into equal four blocks and apply it over and over 
again sequentially since the resulting dynamics is governed by deterministic multiplicative process instead. 



IV. SCALE-FREE PROPERTIES OF WPSL 



Defining each step of the algorithm as one time unit and imposing periodic boundary condition in the simulation 
we can immediately write an exact expression for the total number of blocks as a function of time: N(t) = 1 + 3t. 
However, unlike the square lattice where every block shares its borders exactly with four neighbours, the blocks in the 
WPSL can have their common borders with at least four or more blocks albeit all the blocks have exactly four sides. 
In fact, we can characterize individual blocks of the WPSL by their coordination numbers, which is defined as the 
number of neighbours with whom they have common borders, which we find neither a constant like the square lattice 
nor has a typical mean value like the Voronoi lattice. Instead, the coordination number that each block assumes 
in the WPSL is random and evolves with time. The coordination number disorder in the WPSL can, therefore, be 
regarded as of annealed type. Our data extracted from numerical simulation suggest that the number of blocks Nk(t) 
which have coordination number k (or k neighbours) continue to grow linearly with time Nk(t) = uikt (see figure 
7). On the other hand, the number of total blocks N(t) = J2 k Nk(t) at time t in the lattice also grow linearly with 
time N(t) = 1 + 3i and hence N(t) ~ 3t in the long-time limit. The ratio of the two quantities Pfc(i) = Nk(i)/N(t) 
that describes the fraction of the total blocks which have coordination number k is Pk(t) = rrik/3 and it is clearly 
independent of time vis-a-vis of the size of the lattice (see figure 8). It implies that we can take a lattice of sufficiently 
large size and study its properties as these properties without being worried about its exact size. 
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FIG. 7: Shown is the growth of the number of blocks N k (t) with exactly k — 4, 5, 6, 7 neighbours as a function of time. 
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FIG. 8: Fraction of the blocks pk{t) with k = 4, 5, 6, 7 neighbours are drawn as a function of time. 



We now assume that that blocks in the WPSL are equivalent regardless of their sizes and then ask: What is the 
probability that a block picked at random with equal a priori probability from all the available blocks of the lattice 
have coordination number fc? Answer to this question lies in collecting data for pt as a function of k and plotting 
a normalized histogram which effectively gives the coordination number distribution function p t {k) = P(k) where 
subscript t indicates fixed time. We first plot P{k) in figure 9 to have a glimpse of its behaviour as a function of fc. A 
closer look into the data reveals that there exist scarce data points near the tail. The long tail with scarce data points 
turns into a noisy or fat-tail when we plot lnP(fc) vs ln(fc). This is drawn in figure 10 where data points represents 



12 



0.5 



0.4 



0.3 



0.2 



0.1 - 




FIG. 9: Shown is the coordination number distribution function P(k) vs k which is also equivalent to the degree distribution 
of the DWPSL. 




FIG. 10: Drawn is plot of \n(P(k)) against ln(fc) for the DWPSL network where data points represent average of 50 independent 
realizations. The line have slope exactly equal to 5.66 revealing power-law degree distribution with exponent 7 = 5.66. 



ensemble average over 50 independent realizations and it clearly suggests that fitting a straight line near the tail is 
possible. It implies that the coordination number distribution decays obeying power-law 



P(k) ~ fc-T, 



(50) 



with heavy or fat-tail reflecting scarce data points near the tail-end of figure 10. However, the noise in the tail-end 
complicate the process of identifying the range over which the power-law holds and hence of estimating the exponent 
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7. One way of reducing the noise at the tail-end is to plot cumulative distribution P(k' > k) which is related to 
degree distribution P(k) via 

m _ (5l) 

We therefore plot ln(P(fc' > k)) vs ln(fc) in figure 11 using the same data of figure 10 and find that the heavy tail 
smoothes out naturally where no data is obscured. The straight line ht of hgure 11 has a slope 7 — 1 = 4.66 which 
indicates that the degree distribution (hgure 10) decays following power-law with exponent 7 = 5.66. A power-law 
distribution function is often regarded as scale-free since it looks the same whatever scale we look at it. 
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FIG. 11: Cumulative degree distribution P(k' > k) is shown using the same data as of figure 4. The dotted line with slope 
equal to 7 — 1 = 4.66 is drawn to guide our eyes. 

We thus find that in the large-size limit the WPSL develops some order in the sense that its annealed coordination 
number disorder is scale-free in character. This is in sharp contrast to the quenched coordination number disorder 
found in the Voronoi lattice where it is almost impossible to find cells which have significantly higher or fewer 
neighbours than the mean value k = 6 [l9j |. In fact, it has been shown that the coordination number distribution of 
the Voronoi lattice is Gaussian in character instead. The power-law coordination number distribution in the WPSL 
is reminiscent of the scale-free degree distribution of complex network. The past decade have witnessed a surge of 
interest in the theory of complex network resulting in the dramatic advances in the field thanks to the seminal work 
of A.-L. Barabasi and his co-workers. It is worth mentioning that the dual of the WPSL, obtained replacing each 
block with a node at its centre and common border between blocks with an edge joining the two nodes (see figure 12 
where we have illustrated the network aspect of the WPSL) is a network whose degree distribution has the same data 
as the coordination number distribution of the WPSL. 

The square lattice, on the other hand, is self-dual and its degree distribution is P(k) = 5(k — 4). Further, it is 
interesting to point out that the exponent 7 = 5.66 is significantly higher than usually found in most real-life network 
which is typically 2 < 7 < 3. This suggests that in addition to the PA rule the network in question has to obey 
some constraints. For instance, nodes in the WPSL are spatially embedded in Euclidean space, links gained by the 
incoming nodes are constrained by the spatial location and the fitting parameter of the nodes. Owing to its unique 
dynamics this was not unexpected. Perhaps, it is noteworthy to mention that the degree distribution of the electric 
power grid, whose nodes like WPSL are also embedded in the spatial position, is shown to exhibit power-law but with 
exponent 7 poW er = 4 [20| . 

The power-law degree distribution P(k) has been found in many seemingly unrelated real life networks. It implies 
that there must exists some common underlying mechanisms for which disparate systems behave in such a remarkably 
similar fashion [21J. Barabasi and Albert argued that the growth and the PA rule are the main essence behind the 
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FIG. 12: Nodes in the circle and their links with other nodes illustrate the network topology at step 5 of Fig. 2. New nodes (17, 
18 and 19) and links to be added and links to be removed (broken lines) in step 6 again of Fig. 2 are also shown to illustrate 
the dynamics. 



emergence of such power-law. Indeed, the DWPSL network too grows with time but in sharp contrast to the BA 
model where network grows by addition of one single node with m edges per unit time, the DWPSL network grows 
by addition of a group of three nodes which are already linked by two edges. It also differs in the way incoming nodes 
establish links with the existing nodes. 

To understand the growth mechanism of the DWPSL network let us look into the jth step of the algorithm. First 
a node, say it is labeled as a p , is picked from the (3j — 2) nodes preferentially with respect to the fitness parameter 
of a node (i.e., according to their respective areas). Secondly, connect the node a p with two new nodes (3j — 1) and 
(3j + 1) in order to establish their links with the existing network. At the same time, at least two or more links of 
a p with other nodes are removed (though the exact number depends on the number of neighbours a p already has) in 
favour of linking them among the three incoming nodes in a self-organized fashion. In the process, the degree k p of 
the node a p will either decrease (may turn into a node with marginally low degree in the case it is highly connected 
node) or at best remain the same but will never increase. It, therefore, may appear that the PA rule is not followed 
here. A closer look into the dynamics, however, reveals otherwise. It is interesting to note that an existing nodes 
during the process gain links only if one of its neighbour is picked not itself. It implies that the higher the links (or 
degree) a node has, the higher its chance of gaining more links since they can be reached in a larger number of ways. 
It essentially embodies the intuitive idea of PA rule. Therefore, the DWPSL network can be seen to follow preferential 
attachment rule but in disguise. 



V. SUMMARY 



To summarize, we proposed a model that generates weighted planar stochastic lattice (WPSL) which emerges 
through evolution. Unlike the regular lattice where every block or cell has exactly the same coordination number, the 
WPSL is seemingly disordered both in terms of the coordination number and in terms of its block size distribution. 
One of our primary goal was to find some order in such seemingly disordered system since finding order is always an 
attractive proposition for physicists. To this end, our first attempt was to quantify the size disorder of the lattice 
and we found several interesting results. One of the most interesting results that we found is that x"~ 1 yf^ n 1 
remains a constant V n regardless of the size of the lattice where blocks are labeled as i — 1,2,..., A and Xi and yi are 
their respective length and width. Yet another interesting fact is that the numerical values of the conserved quantities 
x 7~ 1 yt^ n 1 f° r a gi ven value of n except n — 2 are never the same in different realizations which is clearly an 
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indication of wild fluctuation. Of the conservation laws, we found a special one obtained by choosing either n = 1 or 
n = 3 to give X^Li Vi or ^2f=i x i which can be picked as a measure. It implies that if the blocks are populated with 
a fraction of the measure equal to cubic power of their respective length or width then its distribution on the WPSL 
is multifractal in nature - a property that quantifies the wild fluctuation of block size distribution of the WPSL. That 
is, the probability distribution of cc? is multifractal in the sense that a single exponent is not enough to characterize 
its distribution on the support, namely the WPSL, instead we need a spectrum of exponent f(a). We derived an 
exact expression for the /(a) spectrum and obtained information dimension D\ =6/5 with which the entropy of the 
measure scales S(S) <~ In 5~ Dl . To look for a possible origin of mult ifr act ality we also studied the kinetic square lattice 
which is actually the deterministic counterpart of the WPSL. It led to the conclusion that as soon as randomness is 
ceased from the definition of the generator the resulting lattice can no longer be quantified as multifractal rather it 
possesses all the properties of the square lattice but only in the long-time limit. 

Our second attempt was to quantify the coordination number disorder. We have shown numerically that the 
coordination number disorder is scale- free in character since the degree distribution of its dual (DWPSL), which is 
topologically identical to the network obtained by considering blocks of the WPSL as nodes and the common border 
between blocks as links, exhibits power-law. In other words, the coordination number distribution function of the 
WPSL decays following power-law with exponent 7 = 5.66. However, the novelty of this network is that it grows by 
addition of a group of already linked nodes which then establish links with the existing nodes following PA rule, though 
in disguise, in the sense that existing nodes gain links only if one of their neighbour is picked, not the node itself. 
Finally, such multifractal lattice with scale-free coordination disorder can be of great interest as it has the potential 
to mimic disordered medium on which one can study various physical phenomena like percolation and random walk 
problems etc. One may also study the phase transition and critical behaviour in such multifractal scale-free lattice. 
Indeed, phase transition and critical behaviour in complex network have attracted much of our recent attention. The 
advantage of the WPSL over many known complex networks is that its nodes or blocks are embedded in spatial 
position. That is, one may assume that interacting particles located on the sites with great many different neighbours 
and therefore its result must differ from that of a system where sites are located on a regular lattice. We intend to 
work in these directions in our future endeavour. 

NIP gratefully acknowledges support from the Bose Centre for Advanced Study and Research in Natural Sciences. 
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